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Abstract 

We investigate Ad SU(2) lattice gauge theory with Regge-Einstein quantum 
gravity on a dynamically coupled Regge skeleton. To overview the phase dia- 
gram we perform simulations on a small 2 • 4 3 system. Evidence for an entropy- 
dominated disordered, an entropy-dominated ordered and an ill-defined region 
is presented. 
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In recent years various exploratory simulations of pure 4c? quantum gravity have 
been carried out, see ref. [[I], || for recent overviews. In ref. one of the authors 
reported numerical evidence for an "entropy-dominated" phase, where the partition 
function is well-defined, although the Regge-Einstein action is unbounded. Essen- 
tially, these results were confirmed by Beirl et al.,Q, where also some universality 
under changes of the measure was reported. The question of possible physical rel- 
evance of the entropy-dominated phase was first addressed in ref. ||, and it was 
suggested to show that the entropy-dominated phase can accommodate hadronic 
masses many order of magnitude smaller than the Planck mass mp. SU(2) lattice 
gauge theory was suggested to represent the matter fields, because it stands out as 
the simplest 4d asymptotically free quantum field theory. In ref. ||, see ref. [F| for 
a concise summary, a research program along this line of thought was outlined, and 
supported by exploratory numerical results from an 2 -4 3 lattice. For the coupled sys- 
tem the entropy dominated phase survives, and within it an order-disorder transition 
could be located. In a continuum or almost continuum limit, if it exists, the transi- 
tion will have the usual interpretation of a deconfining phase transition, the ordered 
phase defines a quark-gluon plasma and in the disordered phase quarks are confined. 
To provide modest evidence for a consistent picture, one has to show that m^/m-p 
can be decreased. This will require a finite size scaling analysis, and simulations on 
4 • 8 3 and larger systems. They are out of reach for present medium-sized computer 
systems. In contrast, a study of the phase diagram in the (mf>, (3) plane is feasible. 
To map out the phase diagram we performed new simulations, which are reported in 
this paper. 

Let us recall euclidean quantum gravity on a Regge skeleton. The Regge-Einstein 
action is given by 

Sr E = 2ml^2a t A t . (1) 

t 

Here mp is the Planck mass, the sum is over all triangles (d— 2 simplices) of a 4d Regge 
skeleton, A t denotes the area of the corresponding triangle and a t its deficit angle. 
As in previous work, our Regge skeleton will be given by the standard tesselation 
of the hypercubical 4d lattice. It is well-known that the Regge-Einstein action is 
unbounded. Nevertheless the theory could exist due to entropy effects. For volume 
V = const., numerical evidence was reported in ref. || for an entropy-dominated 
region, such that the partition function 




exists for mp < m^ Ax . Here m^ AX limits the well-defined region. As in ref. || we use 
the expectation value < v p >, the average volume of a 4-simplex, to define our length 
unit as / =< v p > l ^ A . Then m^ AX > 0.02/q 2 is found. Henceforth, all quantities will 
be expressed in units of Iq, which will not be explicitly written down anymore. In a 
fundamental length scenario|| Iq stays finite in physical units (for instance Fermi), 
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whereas Iq for a continuum theory. A fundamental length theory is attractive 
under the assumption that it will define a "natural" regularization scheme, such that 
no renormalization occurs when physical quantities are calculated in this scheme. If 
this is not the case, the continuum scenario is preferred, as it offers the possibility 
that (after renormalization) physical results will not depend on microscopic details of 
the regularization, i. e. be universal in the usual sense of lattice gauge theory. These 
two scenarios are presumably indistinguishable in present day numerical simulations. 
Classically our condition V = const, corresponds to a cosmological constant 



\t4V~ 1 J d 4 x^R. (3) 



A = 

2 

This is seen by applying the variational principle to S = mpV^ 1 ^ 2 J d 4 Xy/gR and 
comparing the result with the one obtained from S = mp J d 4 Xy/gR + XV. Equation 
(3) is consistent with the observed value A < 10~ 122 mp, and this would survive 
quantization if renormalization is not necessary. 
Our complete action is 

(4) 

where 

W = -^E^ Re P>(i-^)]- (5) 

Here U t is the ordered product of SU(2) matrices attached to the links of triangle t. 
Coupling of the gauge and the gravity part is entirely achieved through the dynamical 
weight factors Wt- To obtain the correct classical limit, the weight factors Wt have 
to fulfill the constraint 

J2 W t A t = const. V, (6a) 

t 

which should be supplemented by ref. || 

W t > 0. (66) 

These conditions leave ample freedom for choice, and in the following we use the 
barycentric implementation. At mp = 0.005 this system was investigated for various 
(3 values. This is a data point well inside the entropy-dominated region. In the present 
paper we report additional results at mp = 0.010, 0.015, 0.020, 0.025 and 0.0375. 

Let us first consider mp = 0.010, 0.015 and 0.020. There is evidence for the 
disorder-order phase transition at all of these Planck masses and Polyakov loop his- 
tograms for some relevant /3 values are shown in figures 1, 2 and 3. As in ref. || 20,000 
thermalization sweeps were performed before measurement. The subsequent statis- 
tics for each histogram is (3 dependent and lies in the range 20,000-60,000 sweeps. 



More complete numerical results and additional data points are given in ref. |10 
A three-peak structure is sometimes observed when both confined and deconfined 
phases coexist, for example in figure 1. This is suggestive of a first order transition in 
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Figure 1: Histograms of Polyakov loops at two different gauge couplings for Planck 
mass m 2 P = 0.010. 



the vicinity of the coupling. To substantiate this conjecture, would require finite size 
scaling investigations (see for instance ref. [ II]) on 2 ■ 6 3 , 2 • 8 3 , ... systems, and has 



to be postponed until more powerful computational facilities are readily available. 
As function of mp the region of j3 around which the transition occurs does not vary 
much. Tentatively, we identify it to be between (3 = 1.46 and (3 = 1.57. 

The next step is to proceed towards higher Planck masses. At a high enough 
value of the Planck mass one does expect that the system will find itself in unstable 
states, similar as in the pure gravity caseg, |J. Indeed, suggestions to that effect 
begin to appear at m 2 P = 0.025 and (3 = 1.42. Figure 4 shows the time history 
of the gravity action Sre- A two-peak structure, seen usually in first-order phase 
transitions, develops. However, it is unclear whether one can characterize the behavior 
as such because, as in the pure gravity case, the partition function likely ceases to 
exist when one is outside the entropy-dominated region. In any case the two peak 
structure characterizes competing effects of the entropy versus the unboundedness of 
the Regge-Einstein action. 

At Planck mass 0.025 we seem to have reached the region where entropy effects 
from the measure begin to lose control over the unboundedness of the action. In 
figure 6 we present the time history of the gravity action at an even higher Planck 
mass 0.0375 and (3 = 1.42. It is constantly on the rise and the system is clearly 
unstable. 

The (m 2 P , (3) phase diagram is summarized in figure 6. Regions A and C constitute 
the entropy-dominated region which is estimated to be < mp < 0.02, possibly up 
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Figure 2: Histograms of Polyakov loops at two different gauge couplings for Planck 
mass m 2 P = O.Of 5. 
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Figure 3: Histograms of Polyakov loops at two different gauge couplings for Planck 
mass m 2 P = 0.020. 
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Figure 4: Time history of the gravity action Sre for mp = 0.025 and (5 = 1.42. 
Each point is an average over ten consecutive measurements. 



Figure 5: Time history of the gravity action Sre for mp = 0.0375 and (3 = 1.42. 
Each point is an average over ten consecutive measurements. 
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Figure 6: The phase diagram of the system, obtained from simulations on a 2 • 4 3 
system. 

to mp < 0.03. In region B the partition function ceases to exist and is assumed to 
be meaningless. Region C is the ordered phase where, in the infinite volume limit, 
Polyakov loops assume non-zero expectation values and region A is the disordered 
phase where the expectation value of the Polyakov loops are zero. In the present 
paper the hadronic mass (in our unit Iq 1 ) is given by the "deconfining" temperature 

T c = - 1 - , (7) 
2 < xi > 

where < x\ >o is the expectation value of the link length in the short (L = 2) 
direction of the hypercubic lattice. Fairly independent of mp this value stays around 
< xi >o ~ 3. JTTJ Assuming m^ AX = 0.025, one obtains then m-p/m^ < 0.95 for the 



2 ■ 4 3 lattice. With increased computational power one may first envision a finite size 
scaling analysis of 2-L 3 systems towards L — > oo, to confirm and significantly improve 
the accuracy of the phase diagram of figure 6. Simultaneously, one should study an 
4 ■ 8 3 lattice on the present heuristic level, to obtain evidence that the ratio mp / 
can indeed be increased (naively by a factor of two). 
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